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ABSTRACT 

Context. We present the XMM-Newton temperature profiles of 12 bright (L x > 4 x 10 44 erg s ) clusters of galaxies at 0.4 < z < 0.9, 
having an average temperature in the range 5 S kT S 1 1 keV. 

Aims. The main goal of this paper is to study for the first time the temperature profiles of a sample of high redshift clusters, to 
investigate their properties, and to define a Universal law to describe the temperature radial profiles in galaxy clusters as a function of 
both cosmic time and their state of relaxation. 

Methods. We performed a spatially resolved spectral analysis, using Cash statistics, in order to measure the temperature in the 
intracluster medium at different radii. 

Results. We extracted temperature profiles for the clusters in our sample, finding that all the profiles are declining towards larger radii. 
The normalized temperature profiles (normalized by the mean temperature T 500 ) are found to be generally self-similar. The sample 
was subdivided in 5 cool-core (CC) and 7 non cool-core (NCC) clusters, by introducing a pseudo-entropy ratio cr = {T tN /Tout) x 
(EMifi I EMoutT 1 an d defining the objects with cr < 0.6 as CC clusters and those with cr > 0.6 as NCC clusters. The profiles of 
CC and NCC clusters differ mainly in the central regions, with the latters exhibiting a marginally flatter central profile. A significant 
dependence of the temperature profiles on the pseudo-entropy ratio cr is detected by fitting a function of both r and cr, showing an 
indication that the outer part of the profiles becomes steeper for higher values of cr (i.e. transitioning towards the NCC clusters). 
No significant evidence of redshift evolution could be found within the redshift range sampled by our clusters (0.4 < z < 0.9). A 
comparison of our high-z sample with intermediate clusters at 0.1 < z < 0.3, showed how both the CC and NCC clusters temperature 
profiles have experienced some sort of evolution. This can be due by the fact that higher z clusters are at less advanced stage of their 
formation and did not have enough time to create a relaxed structure, characterized by a central temperature dip in CC clusters and by 
flatter profiles in NCC clusters. 

Conclusions. This is the first time that a systematic study of the temperature profiles of galaxy clusters at z > 0.4 has been attempted, 
as we were able to define the closest possible relation to a Universal law for the temperature profiles of galaxy clusters at 0. 1 < z < 0.9, 
showing a dependence on both the state of relaxation of the clusters and the redshift. 
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1. Introduction 

Clusters of galaxies are the largest virialized structures in the 
Universe, arising from the gravitational collapse of high peaks 
of primordial density perturbations. They represent unique sign- 
posts where the physical properties of the cosmic diffuse baryons 
can be studied in great details an d used to trace the past history of 
cosm i c structure formation (e. g. Peebles 19 93tlColes & Lucchinl 
l!995HPeacoclll999l:lRosati et : alJl2002l IVoitl[2005h . As a result 
of adiabatic compression and shocks generated by supersonic 
motion during shell crossing and virialization, a hot thin gas 
permeating the cluster's gravitational potential well is formed. 
This gas reaches temperatures of several 10 7 K and therefore 
emits mainly via thermal bremsstrahlung and is easily detectable 
in the X-rays. The gas key observable quantities are its density, 
temperature, and metallicity. Assuming hydrostatic equilibrium, 
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the gas temperature and density profiles allow to derive the to- 
tal cluster mass a nd thu s to use galaxy clusters as cosmological 
probes (e.g. IVoiJ [20 05). Temperature and density profiles can 
also be combined to determine the intracluster medium (ICM) 
entropy distribution, which provides valuable information on the 
cluster thermodynamic history and has proven to be a powerful 
tool i n investigating non-gravitational processes (e.g. lPratt et al.l 
2006). Radial temperature profiles of the hot ICM in galaxy clus- 
ters are therefore very important to study the gravitational pro- 
cesses responsible for large-scale structure formation and non- 
gravitational energy input into the ICM. 



Early measurements of cluster tempera ture profiles were 
obtained by ASCA an d Beppo-SAX (e.g. Ilkebe et al l 1 19971: 
Mark evitch et all [19981 lEttori et al.l 120001: iFinoguenov et all 
boOlb iNevalainen et alJfeOOlt iDe Grandi & Molendi l2002t) . In 
particular, using ASCA data Mar kevitch et alJ (1998) obtained 
temperature profiles for a sample of 32 nearby clusters, which 
showed significant declines with radius between r = 0.1 r w > 
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and r = 0.6r V! >. They also found out that in clusters with- 
out obvious mergers, the radial temperature profiles outside the 
cool cores were similar when normalized to the virial radius. 
Consistent results were obtained using the spatially resolved 
spectrosc opic capabilities of BeppoSA X. In their analysis of 21 
clusters, iDe Grandi & Molendil d2002l) found declining temper- 
ature profiles, in good agreement with Markevit ch et al.l fl998) 
outside the cores [r > 0.15 - 0.2r V! >], although they were less 
peaked at the center. They also found that for r > 0.2rigo, where 
the gas can be treated as a polytrope, the polytropic index de- 
rived for cool-core (CC) clusters is significantly flatter than the 
index derived for non cool-core (NCC) clusters, corresponding 
to 1.20+0.06 for CC clusters and to 1.46+0.06 for NCC clusters. 

More recently, Chandra & XMM-Newton allowed to mea- 
sure with better accuracy and improved spatial resolution the 
temperature prof iles of galaxy clusters, especially in their cen- 
tral regions (e.g . 
lPrattetal.1 120071 



Vikhlinin et al.ll200l l2006 t iBaldi et alj |2007t 
Leccardi & Molendi 2008), without the com- 



plications introduced by the Point Spread Function (PSF) of 
ASCA and Bepp oSAX. In particular, using Chandra data, 
IBaldi et alJ (|2007) were able to study the temperature profiles 
of 12 clusters in the redshift range 0.1 < z < 0.3, dividing 
the clusters into CC and NCC clusters in a systematic fashion. 
They found that the profiles in the inner O.lriso were showing 
a positive gradient kT(r) oc r 41 , with fi ~ 0.25 in CC clus- 
ters and fi ~ in NCC clusters. They also found that the 
outer region profiles were significantly steeper in the NCC sys- 
tems than in the CC syste ms. This behavior i s in agreement 
with the recent findings of lArnaud et al.l (HblO) who derived a 
Universal pressure profile of galaxy clusters (including both CC 
and NCC clusters) in the REXCESS local sample (z < 0.2) 
with low dispersion, especially in the external regions. Steeper 
temperature gradients in the profiles of NCC systems than in 
those of CC systems are indeed expected if this pressure pro- 
file, analogous for both cluster categories, is combined with 
the den sity profiles (steeper i n CC and flatter in NCC clusters, 
see e.g.|Cavaliere et al.ll2 009). Probing a similar redshift range, 
iLeccardi & Molendil (120081) analyzed a sample of «50 galaxy 
clusters observed by XMM-Newton to measure their radial tem- 
perature profiles. In agreement with previous results they found 
a decline of the temperature in the 0.2 - 0. 6ri80 range. Incon- 
trast with lDe Grandi & Molendi d2002l) and lBaldi etaD (l2007l) . 
they did not find any evidence of a dependence in the slope of 
the outer regions from the presence or not of a cool core. They 
also did not find evidence of profile evolution with the redshift 
out to z ~ 0.3. 

Cosmological simulations have also tackled the measure of 
temperature prof iles in galaxy clus t ers. However, the ex isting 
simulations (e.g. lLoken et al.N2002l : iRoncarelli et al.l 120061) are 
not able to account for the physics at small radii (i.e. in the clus- 
ter core regions) and are not considering the evolution of the 
profiles through cosmic time. 

In this work we aim at studying for the first time the tem- 
perature profiles in a sample of galaxy clusters at z > 0.4, i.e. in 
a redshift range where no comprehensive study is present in the 
literature. The plan of the paper is as follows. In § 2 we describe 
how we selected the XMM-Newton sample we chose and the 
data reduction procedure we followed. We describe our spectral 
analysis strategy in § 3, with particular focus on the background 
treatment procedure (§3.1) and the correction we applied to take 
into account the PSF of XMM-Newton (§ 3.2). In § 4 we present 
the results obtained in our analysis, presenting the temperature 
profiles of individual clusters in the sample and the classification 
method we used to divide them in CC and NCC (§ 4.1), their 
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Fig. 1. Redshift distribution of the XMM-Newton galaxy cluster 
bright sample presented in this paper. 

profiles by the overdensity radius r^oo and by the average tem- 
perature r5oo, and the properties of self similarity they exhibit 
(§ 4.2). We also present several fits of the normalized profiles of 
the total sample (§ 4.2.1), of the CC and NCC clusters samples 
individually (§ 4.2.2), and of the external regions of the profiles 
(§ 4.2.3) . In § 5, we discuss our re sults and compare our sample 
with the ILecca rdi & Molendil (2008) intermediate redshift clus- 
ter sample (§ 5.1), in order to attempt to define a Universal law 
for temperature profiles for galaxy clusters at 0.1 < z < 0.9 
(§ 5.2). Our conclusions are summarized in § 6. 

We adopt a cosmological model with Ho = 70 km/s/Mpc, 
Q m = 0.3, and £2a = 0.7 throughout the paper. Confidence inter- 
vals are quoted at l<x unless otherwise stated. 

2. Sample selection and data reduction 

In Table [TJ we present the list of galaxy clusters analyzed in 
this paper. The clusters considered are picked from the XMM- 
Newton archive sample presented in lBaldi et al.l (1201 2i hereafter 
BAL12). Here we selected all the clusters with at least 3,000 
MOS net counts in order to obtain temperature radial profiles 
having a sizeable number of bins (from 3 to 9, depending on 
the cluster luminosity) and reasonable errors on kT {o-krlkT < 
15%). Thus, a sample of twelve galaxy clusters was selected, 
with a redshift range of 0.4 % z ^ 0.9 being, however, half of the 
clusters located at z < 0.5 (FigureQ~|i. Since this sample was sim- 
ply taken from the archive without any selection criteria apart 
from the number of counts available and the redshift, we cau- 
tion that our cluster sample is far from being a complete sample 
and the relative proportion of CC and NCC clusters may not be 
representative of the relative fraction of CC and NCC clusters 
present in the Universe, within the redshift range explored. The 
full details of data reduction are described in BALI 2 here we 
summarize only the main steps used in processing the data. 

The Observation Data Files (ODF) were processed to pro- 
duce calibrated event files using the XMM-Newton Science 
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Table 1. XMM-Newton archival observations of the 12 bright galaxy clusters analyzed in this paper. 



Cluster 




Obs. Date 


Obs. ID 


tclean.MOS 1 




1 clean, pn 


Nu 


No. of r&ditil bins 










^KsCL. ) 


^KsCL. } 




CIO 20 rm~ 2 1 


CUI1S1UC1CU 


A851 


0.407 


2000 Nov 06 


0106460101 


41.1 


41.1 


30.3 


1.0 


7 


RXCJ0856. 1+3756 


0.41 1 


2005 Oct 10 


0302581801 


24.2 


23.9 


14.7 


3.2 


4 


RXJ2228.6+2037 


0.412 


2003 Nov 18 


0147890101 


24.6 


24.4 


19.3 


4.3 


9 


RXCJ1206. 2-0848 


0.440 


2007 Dec 09 


0502430401 


29.1 


28.8 


20.6 


3.7 


8 


IRAS09 104+4 109 


0.442 


2003 Apr 27 


0147671001 


12.1 


12.3 


8.3 


1.4 


4 


RXJ1347.5-1145 


0.451 


2002 Jul 31 


0112960101 


32.6 


32.4 


27.5 


4.9 


9 


CLJ0030+2618 


0.500 


2005 Jul 06 


0302581101 


15.0 


14.0 


- 


3.7 


3 






2006 Jul 27 


0402750201 


27.0 


27.4 


19.6 










2006 Dec 19 


0402750601 


28.5 


28.6 


20.6 






MS0015.9+1609 


0.541 


2000 Dec 29 


0111000101 


30.6 


30.1 


20.4 


4.0 


7 






2000 Dec 30 


r»i i 1 nAAi r\ 1 

011 1000201 


5.5 


5.4 








MS0451. 6-0305 


0.550 


2004 Sep 16 


0205670101 


25.5 


26.1 


19.4 


3.9 


5 


MACSJ0647.7+7015 


0.591 


2008 Oct 09 


0551850401 


52.9 


53.7 


32.3 


5.4 


6 






2009 Mar 04 


0551851301 


33.2 


34.6 


18.4 






MACSJ0744.9+3927 


0.698 


2008 Oct 17 


0551850101 


40.4 


41.4 


22.4 


5.7 


6 






2009 Mar 21 


0551851201 


63.8 


67.0 


35.8 






CLJ1226.9+3332 


0.890 


2001 Jun 18 


0070340501 


10.3 


10.7 


5.3 


1.8 


4 






2004 Jun 02 


0200340101 


67.2 


67.5 


53.4 







Analysis System (SAS vl 1 . . 1). As described in BAL12, the 
pn data were not used in the spectral analysis, because of prob- 
lems in the characterizati on of the background in pn spectra of 
extended sources (see e.g. lLeccardi& Molendi 2008, and XMM- 
ESAS cookboolfl) and inconsistencies with the measures of tem- 
perature and abundance with the MOS detectors. 

The soft proton cleaning was performed using a double filter- 
ing process. First, we extracted a light curve in 100s bins in the 
10-12 keV energy band by excluding the central CCD, applied 
a threshold of 0.20 cts s _1 , produced a GTI file and generated 
the filtered event file accordingly, in order to remove the ma- 
jority of the flares. The remaining softer flares (not contributing 
significantly to the emission at E > 10 keV) were removed ex- 
tracting a light curve in the 2-5 keV. The resulting histogram was 
fit with a Gaussian distribution having a mean count rate /j, and 
the standard deviation <x. A 3cr clipping algorithm (excluding 
all the time intervals with a count rate larger than + 3cr) and 
visual inspection of the light curves were then used to remove 
the remaining background flaring periods. In Table[T]we list the 
resulting clean exposure times for the MOS and pn detectors. 

After further filtering was applied to remove events with low 
spectral quality, we merged the event files from each detector to 
create a single event file and create a 0.5-8 keV band image of 
each cluster. The resulting image was used to define the extrac- 
tion regions for the temperature profiles and to create a list of 
point sources (using both the EBOXDETECT task and a visual 
inspection of the field to check for false detections). In the case 
a point source is located inside an extraction region, a circular 
region of radius ^ 15" (depending on the source brightness), 
centered at the position of the point source, is excised from the 
spectral extraction region. As in BAL12, the pn data were used 
to create the merged image in order to detect the point sources 
to remove from the cluster spectra at a lower flux limit. 

3. Spectral analysis strategy 

The emission from each cluster was subdivided in several con- 
centric annular regions, allowing each annulus to have a total of 
at least 1,000 net MOS counts and a minimum thickness of 15". 



ftp://xmm.esac.esa.int/pub/xmm-esas/xmm-esas.pdf 



The number of annuli considered for each cluster is reported in 
TableQ] 

3. 1. Background spectral modeling 

The background in our MOS spectra has been modeled 
(in stead of subtracted) follo wing the procedure developed 
by iLeccardi & Molendil (I2008I) and described extensively in 
BAL12. This procedure is preferrable respect to a direct back- 
ground su btraction because it allows the use of Cash statistics 
(ICashl [l979h and avoid problems due to the vignetting of the 
background spectra which would need to be extracted at larger 
off-axis angle than the cluster spectra. Here we summarize just 
th e main aspects of this proce dure, further details can be found 
in ILeccardi & Molendil (l2008h and in BAL12. 

The background model considered consists mainly of five 
components: a thermal emission from the Galaxy Halo (HALO), 
the cosmic X-ray background (CXB), a residual from the filter- 
ing of quiescent soft protons (QSP), the cosmic ray induced con- 
tinuum (NXB) and the fluorescence emission lines. The parame- 
ters for this model were first estimated (using Cash statistics and 
considering the 0.7-8 keV band) from a spectrum extracted in an 
annulus located between 10' and 12' from the center of the field 
of view for both MO SI and MOS 2. To fit the cluster spectra, the 
background parameters (and their lcr errors) were then rescaled 
by the area where t he cluster spectra were ext racted. Appropriate 
correction factors (ILeccardi & Molendlll2008L dependent on the 
off-axis angle) were considered for the HALO and CXB compo- 
nents and a vignetting factor (corresponding to 1.858 - 0.078r, 
with r the distance from the center of the annulus) was applied 
to the QSP component. All these rescaled values (and their la- 
errors) were put in a XSPEC model having the same background 
components us ed in the fit of the 10 '- 12' annulus plus a thermal 
mekal model (iLiedahl et al.ll 19951) for the emission of the clus- 
ter having the temperature, the abundance and the normalization 
free to vary. The lcr errors on the parameters were used to fix a 
range where the normalizations of the background components 
are allowed to vary. This model was tailored to each cluster and 
used to fit all the spectra in our sample. A joint MOS 1 plus MOS2 
fit (using Cash statistics and therefore requiring simply minimal 
grouping, to avoid spectral bins with no counts) was performed 
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to increase statistics. We were able to do that since there are 
no calibration problems between the two detectors and an in- 
dependent analysis of the spectra from MO SI and M0S2 led to 
consistent results in the measures of both kT and Z. 

3.2. XMM-Newton PSF correction and spectral fitting 

In the present analysis we also considered the effect that the 
point spread function (PSF) of EPIC XMM-Newton (HEW ~ 
15") has on the determination of the temperature profiles in our 
cluster sample. This is a different approach respect to BAL12, 
where we were interested only in the measure of the abundances 
in three different radial bins (r < 0.15r5oo, 0.15r5oo < r < 
0.4r5oo, r > 0.4r5oo), a measure which we found not to be af- 
fected by the PSF. In the pr esent sample we used the procedure 
developed in XMM-ESAS dSnowden et al]l2004l) that consider 
the crosstalk Ancillary Response Files (ARFs) between two re- 
gions to remove the contribution of each region from the other 
region and viceversa. Therefore, for each spectral annulus, we 
computed the crosstalk ARF relative to the contribution of this 
annulus to the emission in all the remaining annuli and vicev- 
ersa. All the annuli were therefore fitted simultaneously in order 
to remove the effect of the non-negligible PSF of EPIC, an effect 
that is more important in the core regions and especially if there 
are strong gradients in the emission and spectral parameters over 

the field (as it is the case of cool core clusters). 

The spectra were analyzed with XSPEC vl2.7.0 dArnaudl 

119961) and fitted w ith a single-temperature mekal model 
dLiedahl et alj 1 1995b with Galactic absorption (wabs model). 
The fits were performed over the energy range 0.7 - 8.0 keV 
for both of the two MOS detectors, which were fitted simultane- 
ously to increase the statistics. In our spectral fits, we left free to 
vary the temperature, the abundance, and normalization. Local 
absorption was fixed to the Galactic neutral hydrogen co lumn 
density, as obtained from radio data (Kalb erla et alj|2005l) . and 
the redshift to the value measured from optical spectroscopy (z 
in Table [TJ. We used Cash statistics applied to the source plus 
backgrouncQ, which requires a minimal grouping of the spectra 
(at least one count per spectral bin). A goodness-of-fit for the 
best fit models was computed by using the XSPEC command 
goodness, as described in BAL12. Therefore, for each spec- 
trum in the sample we simulated 10,000 spectra from the best fit 
model, whose C-statistic values were compared with the value 
obtained from the original spectrum. For each cluster, we found 
that the percentage of simulations with a value of the C-statistic 
less than that for the data is always below 20%, an indication 
that the one temperature model used to determine the tempera- 
ture profiles is always a good fit. 

3.3. Determination ofr 500 

Our analysis strategy, which involves testing for self-similarity 
in the temperature profiles at z > 0.4, requires the determination 
of the overdensity radius rsoo and of the average temperature 
Jsoo relative to that radius, com puted as descr i bed in BALI 2, 
adopting the formula derived by Vikh linin et al.l d2006l) : 

r 500 hE(z) = 0.792 (jj^) Mpc, (1) 

where E(z) = (fl m (l + z) 3 + Q A ) a5 . 

2 http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/manual/ 
XSappendixStatistics.html 




LEC/CC , HEC/NCC 

0.2 ' ' ' ' ' 

0.2 0.4 0.6 8 1 

CT 

Fig. 2. Distribution of the values of the pseudo-entropy ra- 
tio <x computed using the radial regions adopted in this pa- 
per vs 1 _fhe_values of cr computed using the radial range used 
in iLeccardi et al.1 d2010l). for the interm ediate redshift clus- 
ter sample of ILeccardi & Molendil d2008l) . The red, green, and 
blue squares represent the systems classified in the literature 
as CC, uncertain, and NCC , respectively. The purple squares 
are the clusters classified in the literature as 'dubious' CC. 
The dotted ho r izont al lines represents the boundary chosen by 
ILeccardi et al.1 d2010b to separate among LEC, MEC and HEC 
clusters. The dashed thick vertical line represent the boundary 
we have chosen to distinguish between LEC/CC and HEC/NCC 
clusters. It is evident how our classification method between CC 
and NCC is in ver y good agreement bo th with the literature and 
with the method of Lecc ardi et al.1 (l2010l) . 



For each cluster, Tgoo was estimated extract a spectrum in the 
0. 15rsoo < r < 0.6r5oo radial range, assuming a first guess value 
of ^oo = 5 keV in order to compute an initial value of 7-500- 
The values of Tgoo and rgoo were then evaluated iteratively until 
a convergence to a stable value of the temperature is obtained 
(ATsoo < 0.01 keV between two successive iterations). The val- 
ues of Tscm and rsoo for each of the clusters in exam are reported 
in Table |2] From this Table, it is apparent how we are sampling 
mostly high temperature clusters, with 75% of the clusters in the 
sample having a temperature kT > 6 keV. 



4. Results 

4. 1. CC/NCC classification of the clusters 

Although the spatial resolution of XMM and the redshift of the 
clusters in our sample in principle could represent a complica- 
tion and may not allow a trivial method to classify galaxy clus- 
ters, we divided our sample in CC and NCC clusters based on 
a comparison of the entropy measured in their core with the 
entropy measured outside the central regions. In particular, we 
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Fig. 3. Temperature profiles for the twelve galaxy clusters in our sample ordered in redshift from left to right and from top to bottom. 
The notes on individual clusters are reported in AppendixlAl 



considered the regions corresponding to r < 0. 15/500 and to 
0.15r5oo < r < 0.4rsoo , defining a pseudo-entropy ratio 



<r = (T, n /Tout) x (EMw/EMoutT 113 , 



(2) 



where Tjm and Tout are the temperature measured in the r < 
0.15r5oo region and in the 0.15r5oo < r < 0.4r5oo annulus, re- 
spectively, and EMin and EMout are the corresponding emis- 
sion m easures. This ratio is very similar to the parameter de- 
fined in lLeccardi et al J (12010) to classify their low-z cluster sam- 
ple (z < 0.24), with the exception of the radial range, which in 
their case was r < 0.05 r2oo for the IN region (corresponding 
to ~ 0.08r5oo) and 0.05r2oo < r < 0.2r2oo for the OUT region 
(corresponding approximately to 0.08r5oo < r < 0.3r5oo). The 



choice of a different radial range was motivated by the higher 
redhift of the clusters in our sample which would make the 
IN region smaller (~ 10") than the PSF of XMM-Newton, in 
most cases. To test whether this choice may introduce a bias 
in the classification, we compu ted the pseudo-entropy rat io for 
the cluster sample presented in iLeccardi & Moleridil d2008l) us- 
ing both methods, and compared the results in Figure [2] It is 
quite clear how the pseudo-entropy ratios calculated with the 
two different choices of radial range are showing a good corre- 
lation with low dispersion especially at low values of <x. In both 
methods, <x is ex pected to be hig her for NCC clusters and lower 
for CC clusters. ILeccardi et al. (|2010) divided their sample in 



three different classes: low-entropy core (LEC, corresponding to 
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Table 2. Values of Tsoo and rsoo (expressed both in Mpc and in arcsec, as computed in BAL12) for the XMM-Newton galaxy cluster 
bright sample presented in this paper. The effective radius R ou , and the outer radius R ou , of the outermost radial bin are reported as 
well. The values of the pseudo-entropy ratio cr and the classification of each cluster are also shown. 



Cluster 


r 50 o (keV) 


r 50 o (Mpc) 


rsoo (") 


Rout (rsoo) 


Rout (rsoo) 


cr 


Classification 


A851 


5.66 + 0.27 


0.975 


180 


0.83 


0.95 


0.896-f- 


NCC/HEC 


RXCJ0856. 1+3756 


6.23 + 0.57 


1.023 


187 


0.53 


0.64 


0.723!^ 


NCC/HEC 


RXJ2228. 6+2037 


7.60 + 0.42 


1.137 


208 


0.97 


1.20 


n /-no+n ion 
633 

-0.092 


NCC/HEC 


RXCJ1206.2-0848 


n Qj+0.60 
4 Q1+8$ 


1.361 


239 


0.66 


0.79 


0.413+0 o» 

°- 582 -!!l 


CC/LEC 


IRAS09 104+4 109 


0.886 


155 


0.55 


0.84 


CC/LEC 


RXJ1347.5-1145 


11.23 ±0.43 


1.366 


237 


0.73 


0.84 


0.349-«i 


CC/LEC 


CLJ0030+2618 


4 69 +036 

t- u -0.30 


0.836 


137 


0.47 


0.58 


0.766!;| 


NCC/HEC 


MS0015.9+1609 


9 48 +0 - 57 


1.184 


186 


0.74 


0.91 


0.877-;" 


NCC/HEC 


MS045 1.6-0305 


9 10 +0 - 56 


1.153 


180 


0.47 


0.56 


555+ 0101 


CC/LEC 


MACSJ0647.7+7015 


9 30 +0 - 45 

y J -0.37 


1.138 


171 


0.61 


0.70 


700 +00< " 


NCC/HEC 


MACSJ0744.9+3927 


g , 4+0.34 


0.995 


139 


0.75 


0.86 


0.451! ;| 


CC/LEC 


CLJ1226.9+3332 


10 16 + ^ 


0.998 


129 


0.50 


0.62 


0.687™ 


NCC/HEC 



cr < 0.45), medium-entropy core (MEC, 0.45 < cr < 0.64), and 
high-entropy core (HEC, cr > 0.64) clusters. For our purposes, 
we do not need such a detailed classification, hence we decided 
to classify our clusters in only two categories: the LEC/CC, with 
cr < 0.6, and the HEC/NCC, with cr > 0.6. Figure |2] shows how 
the choice of cr = 0.6 as a boundary between CC and NCC clus- 
ters is reasonable and agrees well with the independ ent me thods 
of clu ster classification in CC and NCC presented in Ettor i et al.l 
(12010 1. Our sample would then be divided in 5 CC clusters and 7 
NCC clusters, indicating a prevalence of non-relaxed NCC sys- 
tems which is fully expected in the redshift r ange in exam (see 
e.g. IVikhlinin et"atll2007t ISantos et al.ll2010h . The values of cr 
for each cluster are reported also in Table|2] In Figure[3] the tem- 
perature profiles derived for each cluster in the sample, using the 
analysis procedure described above, are shown. 

4.2. Self-similarity of temperature profiles 

The values of rsoo and of Tgoo, computed in § l3.3l and reported in 
Table|2] are used to normalize the average distance from the cen- 
ter of each annulus and the corresponding observed temperature, 
respectively. The normalized temperature profiles are plotted in 
Figure [4] At a first glance, a self-similarity of the temperature 
profiles of the clusters in our sample can be observed, although 
with non-negligible scatter. To the aim of quantifying the self- 
similarity, we considered a function , derived from the relation 
introduced by Vikh linin et al.l d2006l) to fit the deprojected tem- 
perature profiles of their low redshift CC cluster sample: 

T(r) = (*/ar+& 1 

Tsoo (Jf/a)°*+l (1 +{x/b) 1 y 3 ^ 

where x = r/r^oo and A is a normalization parameter. The pa- 
rameters a and b describe the core radius of the inner and outer 
regions of the profiles, respectively, while ao and fio represent 
the inner and the outer slope of the temperature profiles, respec- 
tively. The additional parameter is introduced to adjust the 
relative normalization of the inner profile with respect to the 
outer profile. Figure [5] shows how each of the fitting param- 
eter affects the shape of the temperature profiles described in 
Equation [3] The best fit values of these parameters found by 
Vikhli nin et a D d2006h were A = 1.35, a = 0.045, b = 0.6, 
ao = 1.9, /Jo = 0.45, and = 0.45. Because of the spatial 
resolution of XMM-Newton, the parameter a is clearly insensi- 
tive to our data and it was theref ore fixed to the value of 0.045, 
defined by IVikhlinin et alj (1200 6) . The other parameters are ex- 
pected to be different by fitting our data since we are considering 




0.2 0.4 0.6 0.8 1 1.2 



r / r 500 

Fig. 4. Normalized temperature profiles for all the clusters in the 
sample. CC clusters are plotted with filled squares, while NCC 
clust ers are shown as ope n squares (see § 14.2. 11 1. The best fit to 
the I Vikhlinin e t al. (2006) relation, freeing the normalization, is 
drawn as a blue solid line. The blue dotted line represents the 
best fit to the relation in Equation [3] (with = 0.45), while the 
blue dashed line is the best fit of Equation|5]for r > 0.15r5oo- 



projected temperature profiles (and not 3-D profiles) of a sample 
not exclusively built of CC cluster s. Anyhow, the best f it val- 
ues of these parameters obtained by Vikhlin in et alJ (|2()06) were 
used when we needed to freeze one parameter in order to obtain 
meaningful constraints on the others. 

For co mpleteness, we a lso co mputed the projected equiva- 
lent of the Vikhlin in et a i] d2006h relation. The best fit param- 
eters of the projected relation are slightly different from those 
derived for the 3-D profiles, with the most significant variation 
observed in the inner core radius a (from 0.045 to 0.03), to which 
our data are not sensitive. We run the fits of the temperature pro- 
files using also the projected parameters when freezing one of 
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Table 3. Fit results using the functions defined in ^ 14.21 obtained separately for the CC and NCC cluster sample, by fixing a = 0.045 
and by freezing alternately each one of the other free parameters, except the normalization A. The parameter ao defines the inner 
slope, £o defines the relative normalization between the inner and outer regions of the profiles, b defines the core radius of the 
external regions, while f3o is defining the outer slope. 





Fit function 


^/d.o.f. 




A 




ft 




b 


A 




Eq.|2l(«o = 1.9) 


56.1/28 


1.17 


+ 0.12 




0.61 ±0.15 


0.33 


±0.37 


0.35 ± 0.40 




Eq.EKft = 0.45) 


55.8/28 


1.30 


±0.05 


1.14 ±0.59 




0.28 


±0.26 


0.35 ± 0.32 


CC 


Eq.fJO* = 0.6) 


57.0/28 


1.13 


±0.17 


1.89 ±2.76 


0.65 ± 0.28 






0.72 ± 0.29 




Eq.SOSb = 0.45) 


55.8/29 


1.88 


±0.07 


0.41 ±0.16 


0° 


0.30 


±0.07 






Eq.|3|(a = 1.9) 


35.2/36 


1.24 


±0.12 




0.94 ± 0.24 


0.28 


±0.17 


0.29 ±0.12 




Eq.EKft = 0.45) 


35.2/36 


1.65 


±0.09 


0.18 ±0.66 




0.26 


±0.19 


0.30 ±0.10 


NCC 


Eq.fJO* = 0.6) 


37.0/36 


1.17 


±0.04 


0° 


0.94 ±0.13 






0.62 ± 0.08 




Eq.|3]08o = 0.45) 


35.4/36 


1.20 


±0.04 


unconstr. 


0.91 ±0.12 


0.44 


±0.05 





Notes. (a) Parameter frozen by the fit 
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Fig. 5. Temperature p rofiles as described by the relation of 
Vikhl inin et al.l d2006l) (solid black line) and by varying some of 
its parameters. Decreasing the inner slope parameter ao (green 
dashed line) has the effect of creating a shallower temperature 
dip in the center. If the relative normalization £o between the 
inner and the outer profile is increased (blue dashed line), the 
profile transitions towards the typical shape of NCC clusters. A 
decrease in the outer slope fio (red dashed line) is flattening the 
external part of the temperature profile. 

the fit parameters became necessary, finding consistent results 
with respect to the fits performed using the corresponding 3-D 
parameters. These findings show that a slight change in the value 
of a frozen parameter does not affect significantly the best fits of 
the others, therefore we can conclude that the results of the fits 
presented in the forthcoming sections are quite robust. 

4.2.1 . Fits of CC and NCC temperature profiles 

The temperature profiles of galaxy clusters are known to be 
generally self-similar, showing a decline beyond ~ 15% of 
their virial radius (e.g. IViMilinin etaDl2005b iBaldi et all [20071 



iPratt etalj|2007t iLeccardi & Molendill2008h . However, an ap- 
parent difference is present in their center, with the relaxed CC 
clusters showing a temperature drop while such drop is not ob- 
served in the NCC clusters, which are generally not relaxed. 
Moreover, in a few studies based on low and intermediate red- 
shift clusters, the NCC clusters were found to present a steeper 
temperatu re decline in the outer regions respect to the CC clus- 
ters (e.g. lBaldi et alJl2007|). althoug h this still remains a debated 
issue (e.g. ILeccardi & Molendill2008l) . 

As shown in § 14.11 we divided our sample in CC and NCC 
based on their pseudo-entropy ratio cr. In particular, we defined 
as CC clusters the objects having cr < 0.6, while the objects 
with cr > 0.6 were classified as NCC clusters. The two different 
categories were analyzed individually in this Section. Figure [6] 
shows the normalized temperature profiles for CC and NCC 
clusters separately. We also computed an emission weighted av- 
erage for both subsamples, which, as expected, shows the tem- 
perature drop typical of CC cluster in the central regions, while 
a flat (or increasing) temperature profile is observed at the center 
of NCC clusters. 

To quantify the differences between the two subsamples, we 
fitted the relation in Equation [3] separately to CC and NCC clus- 
ter normalized temperature profiles (Table[3]l. Running a fit leav- 
ing all the parameters in Equation [3] free to vary (except for 
a which is insensitive to our data, and was therefore fixed to 
0.045), no meaningful constraint could be obtained both on the 
CC and on the NCC cluster temperature profiles. Hence, it was 
necessary to alternately freeze one of the parameters (except for 
the normalization A) to the value found in VikMi riin et al.l ([2006) 
(reported in § 14.21 ). and allow the other parameters to fit the data. 
In the case of CC clusters, the better fits are obtained by freez- 
ing the parameter £o or the parameter fio, however in the latter 
case the value of the free parameter £o was frozen to zero by 
the fit. For the NCC clusters the better fits to Equation [3] are ob- 
tained when we freeze ao or £o- If we compare the fits with the 
parameter frozen, which are showing one of the lowest x 1 
and reasonable constraints on the free parameters in both sam- 
ple, CC and NCC clusters differ mainly for the best-fit of the 
inner slope parameter ao which assumes a positive value for the 
CC sample (ao = 1.14 ± 0.59), to fit the temperature drop in 
the center, and a value consistent with zero for the NCC sample 
(ao = 0.18 + 0.66), in order to fit their flat temperature profile 
in the central regions. No significant difference in the external 
slope j3o could be detected, being fto = 0.35 + 0.32 in CC clus- 
ters and fio = 0.30 + 0.10 in NCC clusters. It is also interesting 
to notice how the parameter £o in NCC clusters is always found 
to be ~ 1 (when not frozen to 0.45). A value of ^o ~ 1 makes the 
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Fig. 6. Left: Normalized temperature profiles for the CC clusters in the sample. Shaded areas show the error-weight ed mean of the 
normalized temperature with its error (red) and rms dispersion (magenta) in 4 radial bins. The best fit to IVikhlinin et al.l {2006) 
relation freeing the normalization is plotted as a blue solid line, while the blue dotted line and the blue dashed line represent the best 
fits to the relation in Equation[3](with £o = 0.45) and to Equation [5] (for r > 0.15r5oo), respectively. Right: Normalized temperature 
profiles for the NCC clusters in the sample. Shaded areas show th e error-weighted mean of the normalized temperature with its error 
(blue) and rms dispersion (cyan) in 6 radial bins. The best fit to IVikhlinin et al.l ([2006) relation (leaving the normalization free to 
vary), to the adaptation of the relation in Equation[3](with = 0.45), and to Equation[5]are plotted as a red solid, dotted and dashed 
line, respectively. 



first term in Equation [3] approximately equaling unity, with the 
effect of removing the temperature dip in the center, as expected 
in the profiles of NCC clusters. In the successive fits of NCC 
clusters presented in this paper we will adopt a value of = 1 ■ 



4.2.2. General fits of the whole cluster sample 

We fit the relation in Equation [3] to the normalized tempera- 
ture profiles of the whole sample, including both CC and NCC 
clusters. Leaving all the fit parameters free to vary (except for 
a which is insensitive to our data, and was therefore fixed to 
0.045), we obtain that three of them (A, ao and b) are un- 
constrained and able to variate throughout the whole parame- 
ter space without introducing any difference in the value of the 
statistics, with the other two being loosely constrained (£o = 
0.19 + 0.21, /Jo = 0.35 + 1.08). The refore , as we did for the sepa- 
rate fits of CC and NCC clusters (5 14.2. U . we alternately freezed 
one of the parameters (except the normalization A) to the value 
found by Vikhl inin et al.l (120061) (reported in 5 14.21 1. and allowed 
the other four parameters to fit the data. The results of the fits are 
presented in Table|4] The better fits have been obtained when we 
freezed ao or £o (with minimal difference in their respective^ 2 ) 
yielding a best fit value of = 0.56 + 0.11 andcro = 1.39 + 0.46, 
respectively. These fits are giving also fully consistent values of 
the normalization A and of the parameters b and fio, describing 
the external region of the profiles. 

We also adopted a more general approach to test the depen- 
dence of the profiles on the pseudo-entropy ratio cr (that we used 
to classify the clusters in our sample as CC or NCC), and/or 



on the cluster redshift, by fitting the data with a function of the 
distance from the center r, the pseudo-entropy ratio cr, and the 
redshift z. Equation [3] was modified accordingly in the form: 

T _ (jc/0.045) ff + f 1 

(x/omsr + 1 ' {\+{xlb) 2 f' 

f = & ■ <P + {cr,{ = 0), 
a = ao ■ (f>-(cr,z), 



(4) 



> ± (cr,z) = 1 + T] ■ (1 +0") ± (■ (1 +Z) ■ 

In the function (f> ± (cr,z), the parameters 77 and ( take into ac- 
count the dependence of the temperature profiles on the pseudo- 
entropy ratio cr and on the redshift z, respectively (Figure |7). 
Due to the e xpected evolution of the t emperature profiles in CC 
clusters (e.g. lEttori & Br ighenti 2008), ft is expected to increase 
at high redshift in CC systems and therefore we considered the 
function <p+ to model their evolution. On the other hand, a is ex- 
pected to decrease with redshift and the function 0_ was consid- 
ered for the evolution of this parameter. Although in principle £ 
is expected to have a dependence on z, no meaningful constraint 
could be obtained for its evolution and therefore we imposed this 
parameter to have a dependence only on cr (by setting £ = in 
the function (f> + ). 

Initially, we tested the dependence of the temperature pro- 
files only on the parameter cr (an indicator of their state of re- 
laxation) by fitting Equation|4]to the profiles and considering no 
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Table 4. Fit results using the functions similar to IVikhlinin et all (120061) . defined in § 14.21 obtained by fixing a = 0.045 and by 
freezing alternately each one of the other free parameters, except the normalization A which is always left free to vary. The meaning 
of the fit parameters is the same as in Table [3] The parameters 77 and £ take into account the dependence of the profiles on the 
pseudo-entropy ratio cr and on the redshift z, respectively. 



Fit function 


X^A.o.i. 


A a ft b /3 


1 f 


Eq.|3|(<*0 = 1.9) 
Eq.H«b = 0.45) 
Eq.EKi = 0.6) 
Eq.0060 = 0.45) 


107.4/68 
107.2/68 
111.1/68 
109.2/68 


1.25 ±0.09 - 0.56 ±0.11 0.25 ±0.13 0.28 ±0.11 
1.35 ±0.04 1.39 ±0.46 - 0.22 ± 0.09 0.28 ±0.10 
1.13 ±0.05 3.14 ±2.47 0.72 ±0.10 - 0.61 ±0.10 
1.17 ±0.10 2.45 ±2.13 0.66 ±0.15 0.44 ± 0.08 




Eq.|4|(fc = 0.4; ft = 0.45; 
f = 0) 

Eq.g](fo = 0.4;ft = 0.45) 


109.7/68 
110.3/67 


1.17 ±0.06 4.78 ±5.62 - - 0.26 ± 0.09 
1.48 ±0.15 0.61 ±0.44 - - 0.64 ± 0.49 


0.33 ±0.18 

-0.22 ±0.20 0.12 ±0.39 







b 




0.0 0.2 0.4 0.6 0.8 1.0 
redshift 

Fig. 7. Expected values for the function 0+ introduced in 
Equation [4] in the z - cr parameter space. The black lines rep- 
resent the loci where a combination of the values of z and cr 
are expected to yield identical values of <f> + . The dashed red line 
represent the boundary between CC and NCC clusters. 



redshift dependence at this stage (i.e. fixing the parameter £ = 0). 
Similarly to the fits of Equation [3] we froze one of the free pa- 
rameters in order to obtain reasonable constraints on the others. 
In particular, we chose to freeze the value of ft) = 0.45, corre- 
sponding to the case with the lowest value of^ 2 . This fit however 
is not able to put any meaningful constraint on any of the free pa- 
rameters. Therefore we decided to fix also the value of the exter- 
nal core radius to b — 0.4, a value consistent with the results of 
the fits in Tableg] The fit yields a X 2 ~ 109.7 with 68 d.o.f. and 
points out to a mild dependence on cr of the temperature profiles, 
as indicated by the best-fit value of 77 = 0.33 + 0.18 (Table 0]). 
Even taking into account the uncertainties in the measure of 77 
and of the other parameters, it is clear how the outer part of the 
profiles becomes steeper for higher values of cr (i.e. transitioning 
towards the NCC clusters). This behavior is in agreement with 
the results obtain ed in lower redshift galaxy cluster samples (e.g. 
Bald i et alj2007|) and w ith the Universal pressure profile derived 
bv lArnaud et alJ (l2010h in galaxy clusters in the REXCESS lo- 



cal sample (z < 0.2), showing a low dispersion, especially in 
the external regions. An analogous behavior of the pressure pro- 
files for both CC and NCC clusters, combined with that of the 
density profiles, which are known to be flatter at large radii in 
NCC clusters than in CC clusters (see e.g. lCavaliere et alj [2009: 
Eckert et al. 201 1), would indeed yield steeper temperature pro- 
files in NCC than in CC systems. 

We also added a redshift dependence of the temperature pro- 
files, by fitting Equation|4]to the data, leaving this time also the 
parameter f free to vary. Fixing again ft = 0.45 and b = 04, 
we obtain ax 1 = 110.3 with 67 degrees of freedom, i.e. no im- 
provement in the fit quality. The best-fit value of £ is found to 
be consistent with zero (f = 0.12 ± 0.39), pointing out to no 
significant evolution with the redshift, within the redshift range 
considered (0.4 < z < 0.9). In this fit, also the parameter 77 is 
consistent with zero, an indication that introducing also a red- 
shift dependence removed the mild dependence of the profiles 
on cr detected in the previous fit. 

We caution however that our sample is not a complete sam- 
ple at z > 0.4 with representative fractions of CC and NCC clus- 
ters (§[2]) and the value of the parameter 77 may change, if we use 
a different cluster sample than the one currently in exam. 

4.2.3. Fits of the external regions of the clusters 

As shown before, the shape of temperature profiles in galaxy 
clusters is not easy to describe by using a simple function, with 
a low number of free parameters, to fit the observed data. Indeed, 
in the previous subsections the relations we considered have at 
least five free parameters, a fact that makes harder to obtain 
meaningful constraints on the best-fit values and complicates the 
interpretation of the resulting best-fit functions. The modeling of 
the temperature profiles can be simplified if we consider only the 
external regions of the clusters (r > O.I57-500) in a radial range 
where the influence of a possible cool core is minimal. To fit the 
profiles in these regions we considered Equation [4] and imposed 
that the parameters a and £ are frozen to zero. Thuse we obtain: 

™ A 1 (5) 

7500 (l + Oc/W 

where the parameters b and /3 have the same meaning as in 
Equation [4] To simplify the fit and having just two free pa- 
rameters we also fixed the value of the core radius b to 0.4, a 
value which is more appropriate to the best fits found previously 
in § 14.2.21 for our s ample, than the value of 0.6 measured by 
Vikhlinin etail d2006l) in their low-redshift CC clusters. Without 
considering any dependence on cr and/or z, the best fits to the 
data of Equation [5] yield a. x 1 — 73.2 (both with 51 d.o.f.) with 
best-fit values of the normalization A = 1.13 + 0.05 and of 
/3 = 0.37 ± 0.06 (Table|5). 
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Table 5. Fit results for the external regions of the temperature profiles (r > 0.15rsoo) using the power-law functions defined in 
§ 14.2. 31 The parameters A and /3q describe the normalization and the slope of the power-law, respectively. Also in this case, 77 and £ 
are introduced to take into account the dependence of the profiles on cr and z, respectively. 



Fit function 


^/d.o.f. 


F-test prob. 


A fa 


n ( 


Equation |5](£> = 0.4, 77 = 0, ( = 0) 
Equation |5](£> = 0.4, ( = 0) ' 
Equation |5](£> = 0.4) 


73.2/51 
68.3/50 
67.6/49 


96% 
91% 


1.13 ±0.05 0.37 ±0.06 

1.14 ±0.05 0.95 ±0.32 
1.14 ±0.05 0.64 ±0.54 


-0.37 ± 0.08 

-0.53 ±0.33 0.30 ±0.64 



Table 6. Fit results for the external regions of the temperature profiles (r > 0.15r5oo) using the functions defined in § 14. 2. 31 obtained 
separately for the CC and NCC cluster samples. The fit parameters have the same meaning as in Table[5] Again, the dependence on 
cr (expressed by 77) is not considered here since it was already used subdivide the sample between CC and NCC. 







,^/d.o.f. 


F-test prob. 


A fio 


c 


CC 


Equation |5](£> = 0.4, 77 = 0, f = 0) 
Equation |5](/;> = 0.4, 77 = 0) 


42.2/21 
38.5/20 


82% 


1.11 ±0.08 0.39 ±0.12 

1.12 ±0.08 1.66 ±0.92 


-0.48 ± 0.09 


NCC 


Equation |5](£> = 0.4, 77 = 0, { = 0) 
Equation|5](£> = 0.4, 77 = 0) 


28.2/28 
24.4/28 




1.17 ±0.06 0.37 ±0.06 
1.19 ±0.05 0.06 ±0.01 


4° 



Notes. (<I) Parameter frozen to the upper limit allowed by the fit 



Introducing a dependence from cr (i.e. leaving f = in 
Equation [5] and freeing the parameter 77), a mild dependence 
on this parameter could be found, being the value of 77 = 
-0.37 + 0.08). Freeing also the parameter (, no redshift depen- 
dence could be found (£ = 0.30 + 0.64) within the redshift range 
in exam. 

By fitting separately the CC and NCC samples to Equation|5] 
we obtained a x 2 /d.o.f. = 42.2/21 and a x 2 /d.o.f. = 28.2/28, 
respectively, indicating that the function is a good fit for NCC 
cluster profiles and pointing out to more scattered profiles in CC 
clusters. Both profiles are showing consistent values of the ex- 
ternal slope j3 O80 = 0.37 + 0.06 for NCC, /3 = 0.39 + 0.12 for 
CC clusters). 

We tested also the dependence of the profiles of CC and NCC 
clusters on the redshift, considering the relation in Equation [5] 
and fixing 77 = (Table |6j. The fit yielded an improvement in 
the fit quality of CC clusters (F-test probability ~ 82%), while 
no F-test could be computed in the case of NCC clusters because 
the parameter f was frozen to the upper limit allowed by the fit. 
In CC clusters, we detect a significant dependence of the slope 
from the redshift, as shown by the value of ( — -0.48 + 0.09 
significantly different from zero. However the error on the slope 
is quite large (J3o = 1.66 + 0.92), preventing any quantitative 
conclusion on the evolution of the external temperature profiles 
in CC clusters, within the redshift range probed by our sample. 

5. Discussion 

In § |U we presented the normalized temperature profiles of 
our galaxy cluster sample, showing that they are generally self- 
similar at large radii (r > 0.15r5oo). The profiles exhibit a 
clear decline of the temperature with the radius, although a 
non-negligible scatter is present, especially among the CC clus- 
ters. The func t ion w e used to fit the profiles (derived from 
Vikh linin et al.l (120061) ) represents a good fit to the data in most 
cases, with the largest x 2 observed for the CC sample (where 
the scatter is larger). Although, the profiles showed some depen- 
dence from the pseudo-entropy ratio cr (that we used to clas- 
sify the clusters based on their state of relaxation), dividing 
the sample in CC and NCC cluster, no difference in the exter- 
nal slope between CC and NCC clusters was found, in contrast 
with some of the previous works performed a t lower redshift 
dDe Grandi & Molendill200l iBaldi et al.H"2007l e.g). Moreover, 



an attempt to introduce a dependence also on the redshift gave 
inconclusive results. These facts may be mostly due to the sam- 
ple size, made of 12 galaxy clusters (5 CC and 7 NCC clusters), 
for a total of ~70 data points (~30 and ~45 data points for the CC 
and NCC sample, respectively). Also a reduced redshift range 
(0.4 < z < 0.9) may factor in a lack of evidence of temperature 
profile evolution within our sample. 

Both of these issues could be tackled considering also a 
larger sample of clusters, spanning a different redshift range. 
To this aim, the most suitable sample i s clearly the sample 
of 44 XMM-Newton clusters presented in iLeccardi & MolendU 
(2008), which spans a lower redshift range than our sample 
(0.1 < z < 0.3) and whose data were reduced and analyzed using 
a procedure similar to the one adopted in this paper, especially 
concerning the XMM background treatment. 

5. 1. Comparison with intermediate redshift clusters 

The normalized temperat ure profiles for the sample of 
ILeccardi & Mol e ndl ( 2008I) were obtained using the values of 
rsnn reported in Ettori et all ( 120101) and the values of T500 were 
computed using emission weighted averages of the temperatures 
measured in the radial bins comprised in the 0.15-0.6 r^oo ra- 
dial range. We computed the emission-weighted average tem- 
perature profiles of both CC and NCC clusters, in different ra- 
dial bins in order to have a direct comparison between the two 
samples (probing dif f erent r edshift ranges). Since the size of the 
ILeccardi & Molendj] d2008l) sample is quite large (44 clusters), 
we divided it in two subsamples, considering two redshift bins: 
0. 1 < z < 0.2 and 0.2 < z < 0.3. The 0.4 < z < 0.9 redshift range 
is covered by the sample presented in this work. The emission- 
weighted average temperature profiles in the three redshift bins 
considered are shown in Figure [8] 

The profiles of CC clusters are consistent within the errors at 
every radius except for the central radial bin where the high red- 
shift cluster sample shows a shallower temperature dip respect 
to the lower redshift clusters. This fact may indicate a difference 
between CC clusters at z > 0.4 respect to CC clusters located 
at z < 0.3. A possible reason for this is that the formers are at 
less advanced stage of their formation and would not have had 
enough time to accrete enough mass in order to create a deeper 
potential well and therefore establish a more efficient cooling 
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Fig. 8. Top: Average temperature profiles measured for the CC 
clusters presented in this paper and for the CC clusters of 
ILeccardi & Molendil ((2008), examined in three redshift bins: 
0.1 < z < 0.2 (red triangles), 0.2 < z < 0.3 (green squares), 
and 0.4 < z < 0.9 (blue dots). The profiles are consistent to 
be similar except for the central radial bin where the tempera- 
ture dip in the high redshift clusters is less pronounced. High 
redshift clusters are indeed expected to be at a less advanced 
stage of their formation and not grown enough mass in order 
to create a deep potential well and establish an efficient cooling 
in the central regions. We caution however that another concur- 
ring factor may be the difficulty to fully resolve the core at high 
redshift. Bottom: Average temperature profiles measured for the 
N CC clusters presented in thi s paper and for the NCC clusters 
of ILeccardi & Molendil (2008), examined in three redshift bins: 
0.1 < z < 0.2 (red triangles), 0.2 < z < 0.3 (green squares), 
and 0.4 < z < 0.9 (blue dots). Although the average profiles 
are consistent at every radius, the slope of the high redshift sam- 
ple seems slightly steeper than that of the lower redshift samples, 
which is fully expected because of the less advanced stage of for- 
mation (and consequently less massive hot gas halos) of higher 
redshift clusters. 



in the central regions. However, in this case we have to be cau- 
tious on any claim about evolution for two reasons. First of all, 
it may be very difficult to fully resolve the cluster core regions 
at z > 0.4 because of the XMM-Newton PSF. Although a PSF 
correction was applied (§ I3.21 i. the size of the region correspond- 
ing to the central bin (r < O.lrsoo) in some cases could be of 
the order of the XMM-Newton PSF itself. Moreover, selection 
effects could affect this result. Indeed, our high redshift sample 
is made of only 5 clusters and cannot be considered complete in 
any sense (§|2). 

In the case of NCC clusters, the average profiles in the three 
redshift bins are consistent within the errors at every radius. It is 
worth noticing however that the slope of our high redshift sam- 
ple seems to be slightly steeper than that of the lower redshift 
samples, although no quantitative claim can be obtained from 
this plot. In this case the PSF of XMM is not a biasing factor, 



however selection effects may play a role also in NCC clusters 
(§0. 

A more quantitative comparison could be performed by con- 
sidering the relations introduced in § 14.21 to fit the data points 
from the samples in the three different redshift bins explored. 
For the CC clusters, where the main differences are observed in 
the central regions, the most suitable comparison can be made 
using Equation [3] and fixing the value of the relative normaliza- 
tion between the internal and external region of the profiles to 

= 0.45. Since we are not interested in the external regions of 
the profiles, at this stage we can also fix the external core radius 
to b — 0.4. The results of the fits for the three different red- 
shift bins are shown in Table Q The inner slope ao is flatter in 
the high-z sample presented in this paper wi th respect to both of 
the low er redshift samples extracted from ILeccardi & Molendil 
(120081) . although the difference is statistically significant only in 
the comparison with the 0.1 < z < 0.2 redshift bin (> 2cr). It 
is worth noticing however that also this comparison could be bi- 
ased by the XMM PSF and by possible selection effects. There 
is also a significant difference in the external slope Bo, however 
no quantitative claim could be made about this point because 
fixing the external core radius to b = 0.4 may have biased the 
determination of the slope, especially at lower redshift. 

In the NCC clusters the average profiles are showing a 
slightly steeper slope in our high-z sample with respect to the 
intermediate cluster sample of ILeccardi & Molendil (12008). The 
difference could be quantified fitting the data points for both 
samples using a simplified function as described in § 14.2.31 
(Equation [5} and considering only the external regions of the 
clusters (r > 0.15r5oo). We used this function also to compare 
the external slopes at different redshifts in CC clusters, a task not 
attainable with the more complex function in Equation|3] The re- 
sults of the power-law fits are shown in Table[8] In the case of CC 
clusters there is no significant difference in the power-law slope 
between the two intermediate redshift bins, where the value of 
B is consistent within the lcr errors. A difference just above lcr 
is observed instead between the slope of the high-z CC clusters 
and the slope of the intermediate redshift CC clusters. Also in 
the case of NCC no difference in the slope B is observed between 
the two intermediate redshift bins. The value of B in our high-z 
NCC sample is in th is case significantly h i gher t han in the NCC 
clusters at z < 0.3 in ILeccardi & Molendil (12008) sample (~2<x). 
As stated before, this difference could be due to selection effects 
introduced with our sample selection. However, it is very likely 
that it is instead related with the different stages of the evolution 
of the clusters at z > 0.4 respect to clusters at z < 0.3. It is ex- 
pected indeed that the former would be younger clusters with a 
smaller size hot gas halo with respect to the latter, and therefore 
would show lower temperatures at large radii. 

5.2. A Universal law for temperature profiles 



The data points of the clusters in ILeccardi & Molendil ([2008) 
sample can be used together with our sample to help us define 
a Universal behavior of temperature profiles in galaxy clusters 
at 0.1 < z < 0.9 as a function of both cosmic time and state 
of relaxation (as defined by the pseudo-entropy ratio <x). The 
relations introduced in § 14.21 could be used indeed to fit jointly 
the data points from both our high-z cluster sample and t he in- 
termediate redshift sample of Leccardi & Molendi (2008). The 
results are presented in Table [9] For simplicity, we considered 
only the case of Equation [3] where the value of the external core 
radius has been set at b — 0.4. In all cases, the fit quality is 
not good (x 2 red > 4), indicating a consistent scatter in the nor- 
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Table 7. Fit results for the CC clusters in our sample and in lLeccardi & Molendil (|2008) sample considering three different redshift 
bins and using the relation in Equation [3] obtained fixing the value of the relative normalization between the internal and external 
region of the profiles to £o = 0.45 and the values of the external core radius to b — 0.4. The meaning of the fit parameters is the 
same as in Table [4] 



Redshift range 


X^d.o.f. 


A a Q /3 


0.1 <z<0.2 
0.2 < z < 0.3 
0.4 < z < 0.9 


561.7/94 
186.2/51 
56.2/29 


1.14 + 0.01 2.33 ±0.25 0.31 ±0.05 
1 .02 ± 0.02 1 .96 ± 0.40 0.09 ± 0.05 
1.29 ±0.04 1.01 ±0.39 0.51 ±0.09 



Table 8. Fit results for the external regions of the temperature profiles (r > 0.15rsoo) using the function of Equation [5] (defined in 
§ 14.2.3b obtained separately for the CC and NCC cluster samples in three different redshift bins. The fit parameters have the same 
meaning as in Table 





Redshift range 


^/d.o.f. 


A Po 


CC 


0.1 < z < 0.2 
0.2 < z < 0.3 
0.4 < z < 0.9 


61.4/66 
32.4/39 
42.2/21 


1.09 ±0.01 0.24 ±0.03 
1.07 ±0.02 0.18 ±0.03 
1.11 ±0.08 0.39 ±0.12 


NCC 


0.1 <z<0.2 
0.2 < z < 0.3 
0.4 < z < 0.9 


52.4/57 
50.5/52 
28.2/28 


1.09 ±0.02 0.20 ±0.03 
1.08 ±0.03 0.16 ±0.04 
1.17 ±0.06 0.37 ±0.06 



Table 9. Fit results using the functions similar to Vik hlinin et al. (2006) (defined in § 14. 2b for our sample and the lLeccardi & Molendil 
(2008) sample together, obtained by fixing a = 0.045 and the relative normalization xio - 0.45. The meaning of the fit parameters 
is the same as in Table [4] 



Fit function 


X 2 /d.o.f. 


A a £ Po 


1 <T 


Eq.|3](fe = 0.4) 


1719.0/384 


1.15 ±0.01 1.66 ±0.35 0.41 ±0.06 0.49 ± 0.08 




Eq.|4|(/? = 0.4; ( = 0) 
Eq.|4](/7 = 0.4) 


1663.9/383 
1662.6/382 


1.13 ±0.04 3.13 ± 1.04 0.32 ±0.05 0.17 ±0.04 

1.14 ±0.05 3.85 ± 1.79 0.33 ± 0.06 0.16 ±0.04 


0.29 ±0.11 

0.23 ±0.11 0.17 ±0.09 



malized temperature profiles, Introducing a dependence on cr 
(Equation|4j with £ = 0) the fit is clearly improving (F-test prob- 
ability ~ 99.96%). The best fit values indicate that the external 
part of the profiles is steeper for increasing values of cr, i.e. tran- 
sitioning from CC to NCC clusters. If we take into account also 
a dependence on the redshift of the temperature profiles (by fit- 
ting Equation [4] leaving both 77 and £ free to vary) we obtain 
the following relation that can be considered as closest as pos- 
sible to a Universal law for temperature profiles of clusters at 
0.1 < z < 0.9, attainable with the present data: 

(r,cr,z) = 1.14+0.05 • 7— -, 

T 500 (jt/0.045)" + 1 (1 + (x/0A) 2 f 



(6) 



£ = 0.33 + 0.06 -0 + (cr,<r = O), 
a = 3.85 ± 1.79 ■ <t>-{<r,z), 
J S = O.16±O.O4-0 + (o-,z), 
{(f> ± ((r,z)=l +0.23+0.11 •(! + cr) + 0.17+0.09 • (1 + z) . 



This relation shows that the steepening of the profiles with 
increasing cr is still robust, and the behavior with the redshift is 
in agreement with the results of § 15.11 i.e. higher z corresponds 
to steeper profiles. 



6. Conclusions 

In the present work, we analyzed an XMM-Newton sample of 12 
bright(Lx > 4xl0 44 erg s~') galaxy clusters in the redshift range 
0.4 < z < 0.9, having an average temperature kT > 4.5 keV. This 
sample was extracted from the XMM-Newton sample analyzed 



in BAL12, by selecting all the clusters with at least 3000 MOS 
net counts in order to obtain radial temperature profiles with a 
sizeable number of radial bins and reasonable errors on the tem- 
perature (cr kT /kT < 15%). 

Taking advantage of EPIC XMM-Newton high throughput 
and effective area, which makes it an ideal instrument to per- 
form a spatially resolved spectral analysis, this paper aimed at 
a systematic study of the temperature profiles in galaxy clusters 
at z > 0.4, which is not currently present in the literature. The 
results we obtained can be summarized as follows: 

- we extracted temperature profiles for the 12 clusters in 
our sample, with the cluster extension we probed ranging 
between ~ 500 kpc and ~ 1.3 Mpc from the center. All the 
profiles were found to be declining towards larger radii. 

The temperature profiles of the galaxy clusters in our sam- 
ple, normalized by the mean temperature T500, were found 
to be generally self-similar and could be w ell described by 
a func tion obtained adapting the relation of I Vikhlinin et alj 
(2006) derived for lower redshift clusters. 

We divided the sample in 5 CC and 7 NCC clusters, 
by introducing a pseudo-entropy ratio cr and defining a 
threshold cr = 0.6 between CC and NCC clusters with the 
latters having cr > 0.6. The profiles of the two subsamples 
were found to be different mainly in the inner regions, with 
the inner slope parameter assuming a positive value in the 
CC clusters {a = 1.14 + 0.59), to fit the temperature drop 
in the center, and a value consistent with zero in the NCC 
clusters (a = 0.18 + 0.66), to fit their flat central profile. 
The large errors on the measure of the external slope fio 
gave inconclusive results on the differences between the 
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samples. Also fitting the external regions (r > 0.15rsoo) with 
a simpler function yielded no significant difference between 
the slopes of CC and NCC clusters (J3 a = 0.39 + 0.12 and 
jSo = 0.37 + 0.06, respectively). The lack of any significant 
difference between the two samples could be attributed to 
the small sample size, and therefore to the small number of 
data points available for the fit, especially in the case of CC 
clusters. 

- We introduced a function of both r and cr to fit the data 
points of CC and NCC clusters together. In this case, the im- 
proved statistics allowed to detect a significant dependence 
of the temperature profiles on the pseudo-entropy ratio <x, 
showing an indication that the outer part of the profiles 
becomes steeper for higher values of cr (i.e. transitioning 
towards the NCC clusters). This behavior would agree with 
the re sults obtained in l ower redshift galaxy cluster samples 
(e.g. iBaldi et all [2007b and wi t h the Universal pressure 
profile derived by Arnaud et al.1 (1201 Oh in galaxy clusters 
in the REXCESS local sample (z < 0.2), showing a low 
dispersion, especially in the external regions. 

- In all the fits attempted, no evidence of redshift evolution 
could be found within the redshift range sampled by our 
clusters (0.4 < z < 0.9). 

- We compa red our sample with the intermediate cluster 
sample of iLeccardi & Molendil (120081) at 0.1 < z < 0.3, 
finding significant differences in both the CC and the NCC 
cluster samples. In particular, we found that CC clusters at 
z > 0.4 shows a temperature dip in the center less deep with 
respect to CC clusters at z < 0.3, while NCC clusters at 
z > 0.4 showed steeper temperature profiles with respect to 
NCC clusters at z < 0.3. Both these facts can be due by the 
circumstance that higher z clusters are at less advanced stage 
of their formation and did not have enough time to create a 
large massive hot gas halo comparable in size and mass with 
that of lower redshift clusters. 

- We defined for the first time the closest possible relation to a 
Universal law for the temperature profiles of galaxy clusters 
at 0.1 < z < 0.9. This relation shows a dependence on both 
the state of relaxation of the clusters and the redshift. 

Although some of the results obtained in this paper can be 
in principle biased by possible selection effects introduced in 
extracting the sample from the XMM archive (§ 0, we stress 
that this is the first time that a systematic study of the temper- 
ature profiles in galaxy clusters at z > 0.4 has been attempted. 
Additional deep XMM-Newton and Chandra observations, and 
most likely a new generation of X-ray observatories, would be 
needed to improve the current knowledge of the temperature dis- 
tribution in the hot gas of galaxy clusters at high redshift. 
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Appendix A: Notes on individual clusters 

A851 - A851, the lowest redshift cluster in the sample 
(z = 0.407), has been observed by EPIC XMM-Newton for 
-56 ks (ObsID: 0106460101), of which ~ 41 ks were useful 
for scientific analysis after removing the periods of high flares. 
Its temperature profile is quite flat in the central regions with a 
temperature of ~ 6 keV, showing a significant decline only in 
the outer radial bin (r > 700 kpc). As indicated by a value of 
cr = 0.896+qj^, this cluster is classified as NCC. 

RXCJ0856.1+3756 - RXCJ0856. 1+3756 (z = 0.411) has been 
observed by XMM for ~ 29 ksec (ObsID: 0302581801), with 
~ 24 ks available for our analysis after cleaning the event 
files from soft proton flares. Although the temperature profile 
consists of only four radial bins, a hint of a temperature drop 
in the central bin (r < 170 kpc) is observed, although its value 
of cr — 0.723^qj^ puts this cluster in the NCC category. The 
temperature profile then declines from kT as 9 keV down to 
kT x 5.5 keV) in the outer two bins (r > 250 kpc). 

RXJ2228.6+2037 - EPIC XMM-Newton observed the cluster 
RXJ2228.6+2037 (z = 0.412) for a total of ~ 26 ksec (ObsID: 
0147890101), almost entirely useful for our analysis. Although 
a hint of a drop in the temperature is observed in the central 
two radial bins (r < 150 kpc), with kT x 8 keV in the central 
region, this cluster is classified as NCC (cr = 0.633+° ^°). The 
temperatures reaches its maximum (kT « 9 keV) at r ~ 200 kpc 
and then declines smoothly down to ~ 4 keV at radii larger than 
~ 850 kpc. 

RXC J 1206.2-0848 - The cluster RXCJ1206.2-0848 (z=0.440) 
has been observed twice by EPIC XMM-Newton (Obslds: 
0302581901, 0502430401), however one of the two observation 



13 



A. Baldi et al.: Self-similarity of temperature profiles in distant galaxy clusters 



(ObsID: 0302581901) was heavily affected by soft proton flares 
and was therefore discarded. The ~ 30 ks of the other observa- 
tion were almost entirely available for scientific analysis. The 
temperature profile can be considered approximately constant 
around kT a 10 keV, although several fluctuations are observed 
in the profile and the temperature shows a drop of ~ 3.5 keV 
between the center and the adjacent radial bin. This cluster is 
classified as CC, presenting a value of cr — 0.413^ ^1 

IRAS091 04+41 09 - IRAS09 104+4 109 (z=0.442) has been 
observed by EPIC XMM-Newton for a total of ~ 14 ks 
(~ 12 ks of clean exposure time, ObsID: 0147671001). The 
temperature shows a hint of a temperature drop in the center 
(AkT ~ 0.5 keV), and decreases from ~6.5 keV to ~4.5 keV 
between the region located at ~ 100 kpc from the center and the 
one located at ~200 kpc, remaining constant around the latter 
value going towards larger radii. This cluster shows a value of 
cr = 0.582+J5 |g and is therefore classified as NCC. 

RXJ1 347.5-1 145 - RXJ1347.5-1 145 (z=0.451), representing 
the brightest cluster in our sample with its X-ray luminosit y 
of 4.2 x 10 45 erg s 1 (0.1-2.4 keV band. lEbeling et al.ll2010h . 
has been observed for ~ 38 ks (~ 32 ksec after cleaning for 
flares) by XMM-Newton (ObsID: 0112960101). This cluster 
shows evidence a very strong cool core with the temperature 
dropping from kT a 13.5 keV to kT a 10 keV in the central bin 
(r < 90 kpc) and is classified as a CC cluster (cr = 0.349+°°^, 
the lowest value among the clusters in our sample). Although 
some irregularities in the profile are present, the temperature 
then declines smoothly going towards the outer regions, reach- 
ing a value around ~ 6 keV at r > 750 kpc. 

CU0030+2618 - The cluster CLJ0030+2618 (z=0.500) 
has been observed three times by XMM-Newton (Obslds: 
030258 1 1 1 , 040275020 1 , 040275060 1 ) for a total of ~ 94 ksec, 
of which ~ 70 ksec were available after removing the high- 
flaring intervals. The number of counts available did not allow 
to divide the emission in more than three radial bins, with 
the temperature profiles being approximately constant around 
~ 5 keV. CLJ0030+2618 is classified as a NCC cluster, as 
indicated by a cr = 0.766+^07. 

MS0015.9+1609 - MS0015.9+1609 (z=0.541) has been ob- 
served twice by XMM (Obslds: 0111000101, 0111000201) 
for a total of ~ 45 ks (~ 36 ks of clean exposure time). The 
temperature profile declines smoothly from kT « 12 keV 
measured at the center down to a minimum of ~ 7 keV observed 
as far as ~ 900 kpc from the center, making of it a clear example 
of a NCC cluster (cr = 0.877+° |^). 

MS045 1.6-0305 - The galaxy cluster MS045 1 .6-0305 (z=0.550) 
has been observed by XMM-Newton (Obsld: 0205670101) for 
a total of ~ 45 ks. The observation was heavily affected by 
soft proton flares and we were able to use only ~ 26 ks of the 
original exposure. The observed temperature shows a drop in 
the central 200 kpc (kT ~ 9 keV) and is then quite constant 
around ~ 11 keV between 200 and 450 kpc. It then declines to 
~ 6 keV at r > 450 kpc. The value of cr = 0.555+^ measured 
for this cluster puts it in the CC category. 

MACSJ0647.7+7015 - EPIC XMM-Newton observed twice 
(ObsIDs: 0551850401, 0551851301) the galaxy cluster 
MACSJ0647.7+7015 (z=0.591) for a total of ~ 145 ks 
(~ 88 keV of clean exposure time). The temperature ob- 



served in this cluster shows a smooth decline from the center 
(kT m 12 keV) towards the outskirts (kT « 6-7 keV at 
r > 450 kpc), making of it another obvious example of NCC 
cluster (cr = 0.700+^). 

MACSJ0744.9+3927 - Two EPIC observations of the cluster 
of galaxy MACSJ0744.9+3927 (z=0.698) are present in the 
XMM-Newton archive (ObsIDs: 0551850101, 0551851201) 
accounting for a total exposure time of ~ 154 ks (~ 106 ks 
after cleaning for high-flaring periods). MACSJ0744.9+3927 
is another examples of strong CC present in our sample (as 
indicated also by cr = 0.451^ ^), showing a temperature drop 
from ~ 10 keV (at ~ 150 kpc from the center) to ~ 7.5 keV in 
the central radial bin (r < 100 kpc). The temperature smoothly 
decline going towards the outskirts and stabilizes around a value 
of ~ 6 keV at r > 500 kpc. 

CLJ1226.9+3332 - CLJ1226.9+3332, the farthest cluster 
in our sample (z=0.890), has been observed twice by EPIC 
XMM-Newton (Obslds: 0070340501, 0200340101) for a total of 
~ 139 ks. Both observation were heavily affected by soft proton 
flares (in particular ~ 75% of ObsID 0070340501 was affected 
by soft protons), leaving just ~ 77 ks of clean exposure time 
available for scientific analysis. The temperature shows a sharp 
decline from the center (kT w 13 keV) towards the outer radial 
bin as far as ~ 500 kpc from the core (kT * 6 keV). This cluster 
is classified as a NCC (cr = 0.687+" !^). 
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